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Abstract 

Recent studies in mammalian hearts show that left ventricular wall thickening is an important 
, mechanism for systolic ejection and that during contraction the cardiac muscle develops significant 

stresses in the muscular cross-fiber direction. We suggested that the collagen network surrounding 
the muscular fibers could account for these mechanical behaviors. To test this hypothesis we develop 
a model for large deformation response of active, incompressible, nonlinear elastic and transversely 
0^ ' isotropic living soft tissue (such as cardiac or arteries tissues) in which we include a coupling effect 

\ between the connective tissue and the muscular fibers. Then, a three-dimensional finite element 

formulation including this internal pseudo-active kinematic constraint is derived. Analytical and 
i^h ' finite element solutions are in a very good agreement. The numerical results show this wall thickening 

C$ , effect with an order of magnitude compatible with the experimental observations. 

Keywords : Constitutive law, Finite element method, Living tissue, Hyperelasticity, Nonlinear partial 
differential equations, Anisotropic material. 
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\Q ' 1 Introduction 

It is known that the transverse shear along myocardial cleavage planes provides a mechanism for a 
normal systolic wall thickening [14]. Indirect evidences indicate that the characteristics of the passive 
extracellular connective tissue in the myocardium is an important determinant of ventricular function 
([15], [23], [7]). An appropriate constitutive law for the myocardium should therefore incorporate the 
most important features of its microstructure. A sound theoretical formulation for material laws of the 
active myocardium is essential for an accurate mechanical analysis of the stresses in the ventricular wall 
during the whole cardiac cycle. The wall stress distribution is one of the main factors governing the 
myocardial energetic [25], the coronary blood flow [3], the cardiac hypertrophy [23], and the fetal heart 
growth [19]. To date we do not have any reliable technique to evaluate the stress in the cardiac muscle, 
therefore, mechanical models are useful in cardiology to assess the functional capacities of the human 
heart. Several numerical models using a finite element (FE) analysis have been performed to simulate 
the left ventricular performance ([13], [29]). The mechanical behavior of the connective tissue is often 
assumed isotropic [20] . This last assumption is not in agreement with the experimental results obtained 
on a sample of active myocardial rabbit tissue. Lin and Yin [15] showed that, during an active equibiaxial 
stretch test, there are significant stresses developed in the cross-fiber direction (more than 40% of those 
in the fiber direction) that cannot be attributed to nonparallel muscle fibers (MF). 

Therefore, the purposes of this paper are to: (i) suggest a realistic pseudo-active kinematic law coupling 
the passive connective tissue to the muscle fibers, which may explain the developed tension in the cross- 
fiber direction observed by Lin and Yin ([15]), (ii) formulate an active three-dimensional material law for a 
nonlinear hyperelastic and incompressible continuum medium, which takes care of these coupling effects, 
(iii) derive the related three-dimensional finite element (FE) formulation, and (iv) test the accuracy of 
the proposed numerical method. 



00 

o 



- 1—1 

X 



"Laboratoire de Mathematiques, Universite de Savoie, 73376 Le Bourget du Lac, France, email:Christian Bourdarias 
@univ-savoie.fr 

tLaboratoire de Mathematiques, Universite de Savoie, 73376 Le Bourget du Lac, France, email: Stephane. GerbiOuniv- 
savoie.fr 

t Universite Joseph Fourier, Laboratoire TIMC-IMAG, Grenoble, France and Institut de l'lngenierie et de l'lnformation 
de Sante (IN3S) Faculte de Medecine de Grenoble Btiment Taillefer (Bureau B57) 38706 La Tronche Cedex, France, email: 
Jacqucs.Ohayon@imag.fr 



2 Microstructure of the cardiac tissue 



2.1 Muscle fiber organization 

Anatomical observations have shown that the cardiac muscle tissue has a highly specialized architecture 
[26]. This structure is composed primarily of cardiac muscle cells, or myocytes, that are 80 to 100 /im 
in length and are roughly cylindrical with cross-sectional dimensions of 10 to 20 /im. These cells are 
arranged in a more or less parallel weave that we idealize as "muscle fibers" (MF). We shall denote the 
local direction of this group of cells by the unit vector f and refer to it also as the local "fiber" direction 
with the understanding that individual continuous MF do not really exist. Experimental measurements 
have shown that the MF direction field defines paths on a nested family or toroidal surfaces of revolution 
in the wall of the heart [26] . These results show a continuously changing orientation f of the MF through 
the wall, circumferential near the midwall and progressively more inclined with respect to the equatorial 
plane when moving toward either the epicardium or the endocardium. 

2.2 The cardiac connective tissue organization 

Myocytes and coronary blood vessels are embedded in a complex extracellular matrix which consists of 
collagen and elastin, mainly. Caulfield and Janicki [2] used the scanning electron microscope (SEM) to 
reveal the basic organization of this connective tissue network. Their studies on the connective tissue of 
mammalian heart muscle give the description of the extracellular structures and their arrangement relative 
to cardiac muscle cells. They described the three following classes of connective tissue organization: (i) 
interconnections between myocytes, (ii) connections between myocytes and capillaries and, (iii) a collagen 
weave surrounding group of myocytes. When viewed by SEM, groups of myocytes can be seen to be 
encompassed by a rather prominent meshwork of fibrillard collagen, and short collagen struts attach the 
myocytes subjacent to this meshwork to it. 




Figure 1: Collagen network surrounding the myocytes. 



3 Constitutive law in continuum mechanics 

3.1 Coupling between muscle fibers and collagen network 

Extrapolations from muscle fiber arrangement to myocardial stress are realistic when also taking account 
the effect of the connective tissue. We believe that a part of that connective tissue, surrounding group of 
myocytes, is responsible for active tension developed in the perpendicular direction of the muscle fibers 
running on the tangential plane of the ventricular wall. 

Based on the previous SEM observations, we proposed a connective tissue organization illustrated on 
Figure 1. We assumed that the myocytes are roughly cylindrical and that groups of myocytes are 
surrounded by inextensible collagen networks. So, during the contraction, the myocytes diameter increases 
and because the collagen network is inextensible, the adjacent muscle cells become closer. Thus the 
pseudo-active kinematic relation between the muscle fiber and cross-fiber extension ratios (noted A f and 



X c f, respectively) is h(Xf, A c /) = with: 

h(X f , Ac/) = 1 - A c/ + (tt - 2)(1 - A7 1/2 )-| (1) 
with Z) = 4a + d where a is the initial myocyte radius and d is the distance between the two cells. 
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Figure 2: Schematic illustration of the internal pseudo- active kinematic constraint induced by the collagen 
network surrounding the myocytes. A) Before contraction. B) After or during contraction. 

3.2 Constitutive law for the myocardium under internal pseudo-active kine- 
matic constraint 

To be consistent with our mathematical formulation, the letter <I? is used for non elastic gradient tensor 
and the letter F is used for elastic gradient tensor. The activation of the muscle fibers changes the 
properties of the material and at the same time contracts the muscle itself. To have a continuous elastic 
description during the activation of the myocardium, we used an approach similar to the one proposed 
by Chadwick [4], Ohayon and Chadwick [20], Taber [27]. From its passive zero-stress state P, the free 
activation of the muscle fibers is modeled by the following two transformations (Fig. 2): the first one 
(from state P to virtual state ^4o) changes the material properties without changing the geometry, and 
the second one (from A$ to A) contracts the muscle without changing the properties of the material. 
Thus, the former is not an clastic deformation and is described by the gradient tensor &pA a = I where 
I is the identity matrix. In that first transformation, only the strain energy function is modified using 
an activation function (3, which may depend on the cardiac cycle time and some ionic concentration 
(calcium for instance). The second transformation is an elastic deformation caused only by the active 
tension delivered by the fibers and takes care of the internal kinematic constraint (Eq.(l)). This last 
transformation is described by the gradient tensor Fa a- Thus the transformation from state P to state 
A is a non elastic transformation (<&pa = &pa ~Fa a), but can be treated mathematically as an elastic 
one because <&pa — Fa a- Finally, external loads are applied to state A deforming the body through 
into C (Fig. 3). 
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Figure 3: Description of the active rheology approach. 

The change of the material properties of the myocardium during the cardiac cycle is described by a 
parameter-dependent strain-energy function per unit volume of state P noted W(Eph)' 



with 



W(E PH ) = -Iph(MVph) - 1) + W*(E PH ) + 5 ah W pseudo (E PH ) 

Z active 



W*(E PH ) - W pas (E PH ) +0Wi ct (E PH ) 



(2) 



(3) 



where Epg is the Green's strain tensor at an arbitrary state H calculated from the zero strain state P 
(the state H could be one of the states Aq, A or C shown in figure 3), pn is the Lagrangian multiplier 
resulting of the incompressibility constraint det 4>pp = 1 (see [8, 16, 12]), .^(Epp) is the determinant of 
the right Cauchy-Green strain tensor Cpp (Cpp = 2Epp + I), W pas represents the contribution of the 
surrounding collagen matrix and of the passive fiber components, W[ ct arise from the active component 
of the embedded muscle fibers, and (5 is an activation parameter equal to zero at end-diastolic state and 
equal to one at end-systolic state (0 < (3 < 1). The scalar 5 ah is equal to one if state H is the state A 
and zero if the two states H and A are distinct. The term W/ ct (Epp) gives the variation of the muscle 
fibers properties during the cardiac cycle. The pseudo-active strain energy function expressed in the last 
term of the right hand side of the Eq.(2) is introduced in order to satisfy the kinematic condition (Eq.(l)) 
and is given by: 



(E PJJ ) = --q H h(E PH ) 



(4) 



The scalar gp introduced in Eq.(4) serves as an additional indeterminate Lagrange multiplier which 
contributes to the pseudo-active stresses at state H in fiber and the cross-fiber directions, and /i(Epp) 
is the function defined in Eq.(l), which may be rewritten as: 



h(E PH ) = 1 - il 12 + (tt - 2) (1 - i L 



-1/4 



»5 



(5) 



fp • C pp • fp in 



where I4 and 7 6 are two strain invariants given by I 4 (Epp ) = fp • C pp • fp and J 6 (Epp ) 
which the fiber and the perpendicular fiber directions (this last one corresponding to the direction of the 
collagen struts) are respectively characterized in state P by the unit vectors fp and fp. In an arbitrary 
deformed state H, the direction of these two unit vectors are noted fp and i' H and are respectively defined 
by: 



fp = 



'pp 



■f F 



'pp 



•fp 



and f' H = 



*pp • n 



f^ 
r p 



-pp 



The tensor Cpp is the right Cauchy-Green strain tensor (Cpp = 2Epp + I = <frp H $>PH )• The 
superscript 'T' is used for the transpose matrix and || • || stands for the euclidian norm. Note that I4, 
and are directly related respectively to the fiber and cross-fiber extension ratios (we have I4 = X"j and 
h = X 2 c j). In our notations A/ is related to the fiber direction fp and A c / to the cross- fiber direction 



A 



i' H (Figure 2). We treat the myocardium as a homogeneous, incompressible, and hyperelastic material 
transversely isotropic with respect to the local muscle fiber direction. 
In this study, the passive strain-energy function is [15] 

W pas (V PH ) = Cf(e<3 -1) (6) 

with Q = CP(I 1 -3) 2 + Cl(I 1 -3)(h-l) + Cl(h-l) 2 (7) 

For the active strain-energy we modified the function found by Lin and Yin [15] by substracting the 
"beating term" Cg: 

W act {EpH) = C?(h - 3)(/ 4 - 1) + C 2 a (/x - 3) 2 + CI(I 4 - l) 2 + C a A (h - 3) (8) 

where Cf , i = 1, • • • , 4 and Cf , i = 1, • • • 4 are material constants and Ji is the first principal strain 
invariant given by Ii(Eph) = tr Cph ■ 

The beating term is defined as the part of the active strain-energy function responsible for the change of 
geometry when the muscle is activated and submited to no external loading. To incorporate the beating 
behavior, the parameter-dependent beating tension (3T^ was applied in the deformed fiber direction. 
In our approach, the active loaded state C of the myocardial tissue is obtained in two steps. In the first 
step and at a given degree of activation (3, we derived and quantified the internal pseudo-active stresses 
by looking the free contraction configuration of the tissue (state A, Figure 3). Then, in a second step 
we applied the loads on the active myocardial tissue under the internal pseudo-active stresses previously 
found. 

Step 1: determination of the free contraction state A- During the cardiac cycle and at a given 
degree of activation (3, the Cauchy stress tensor in state A (noted ta) is given by: 

T A = - PA I + *PA dW *J; EPA) <t> T PA + 0T<- O) fA ® fA + (9) 
Ot,pA 

pscudo dW pseu do (Epi) 

with T a A cUve = *p A (10) 

where the symbol (g> denotes the tensor product. The postulated mechanical coupling law (Eq.(5)) induces, 
during the contraction, a pseudo-active stress tensor: 

These two stress tensor components T{ and are activation-dependent and behave as some internal 
tensions in the fiber and cross-fiber directions of unit vectors and f' A , respectively. These pseudo- active 
tensions are defined by: 

dW peeudo dW pseudo 

T J = 2 acUve || . fp ||2 . ye/ = g acUve || ^ . fj L ||2 (12) 



Step 2: determination of the physiological active loaded state C- These previously found internal 
pseudo-active tensions Tj{ and T A were introduced in the expression of the stress tensor at loaded state 
C. Therefore, at a given degree of activation /?, the Cauchy stress tensor in the physiological state C 
(noted tc) is given by: 

r c = -Pd + *rc dW Q^ p P c C) * T p C + (> (0) + T A ) f c ® f c + T C J f' c f' c (13) 

The suggested constitutive law for the active myocardium (Eqs.(2)-(13)) allows to simulate the left 
ventricle behavior during the whole cardiac cycle. Thus, in this law: (i) the anisotropic behavior is 
incorporated in the expressions of passive, active and pseudo-active strain energy functions by the terms 
Ii and Iq, (ii) the kinematic contraction is accounted for by a beating tension [3T^ in the fiber direction, 
(iii) the change of properties is expressed by the active strain energy term (3W act , and (iv) the coupling 
effect between the collagen network and the MF is accounted for by the two internal pseudo-active tensions 
T A and T A in the fiber and cross fiber directions fc and f' c , respectively. 
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4 Variational formulation and finite element method 



The undcformcd body state P contains a volume V bounded by a closed surface A, and the arbitrary 
deformed body state is, as before, noted H. The corresponding position vectors, in cartesian base unit 
vectors, arc R = Y R e R and r = y r e r , respectively. However, we write the equations with suitable 
curvilinear systems of world coordinates noted Q A in the reference configuration (state P) and 9 a in the 
deformed configuration (state H): see Fig. A in appendix A. In this paper we use the same conventional 
notations (Table 1 in appendix A) for vectors, tensors and coordinates systems than Costa et al. [6, 7], 
where: 

- Capital letters are used for coordinates and indices of tensor components associated to state P, and 
lower case letters are related to state H. 

- G and g are the base vectors in states P and H , respectively, for which parenthetical superscript 
indicates the associated coordinate system (for example = dH/dX 1 = R^ and g^ = 
dr/dx l ^v { f). 

The Lagrangian formulation of the virtual works principle is given by ([6, 16]) 

J ' Pjj r $j a V i(Su a ) dV = J p(b a -j a )5u a dV + (1 -5 AH ) J s.SudA (14) 

V V A 2 

where P H J are the components of the second Piola-Kirchhoff stress tensor at state H, Ph, referred to 
the base tensor G^ <8> Gj°\ = d9 a /dX I are the components of the gradient tensor 3>p# in the 
base tensor g^P ® G^' 7 , Su = 5u a g^ a is an arbitrary admissible displacement vector, Vi(Su a ) — 
dSua/dX 1 — g^\ • g^^Su/s are the components of the covariant differentiation vector Su in the base 
vectors g^ a (i.e. Vi{6u) = V/(,5u a )g (e)a ). The previous differentiation is done with respect to the 
locally orthonormal body coordinates ( X 1 , / = 1,2,3) which coincide with the local muscle fiber 

(8) 

direction in state P. The material density in the undcformcd body state P is p, b = b a g a is the body 
force vector per unit mass, 7 = 7 Q g« is the acceleration vector, s is the surface traction per unit area of 

A, and A2 is the part of A not subject to displacement boundary conditions. The Lagrangian formulation 
for incompressibility is given by 

jf(det fl #-l) P *dV = (15) 

(x) 

where the metric tensor 5} j is defined in table 1, and p* is an arbitrary admissible pressure. Lastly the 
Lagrangian formulation for the additional pseudo-active kinematic constraint is given by 

S AH f h(I 4 ,I 6 )q*dV = (16) 
Jv 

for all admissible q* . Eqs.(14)-(15) -(16) represent the variational formulation of a system of nonlinear 
partial differential equations. For an incompressible medium (det &ph = 1), the relation between the 
second Piola-Kirchoff stress tensor P# and the Cauchy stress tensor th is ([16]) 

P fl =^.r fi .(^f (17) 

A complete expression of the components of Ph in both states H = A and H = C are given in appendix 

B. The surface traction per unit of undcformcd area of A, s = s a ga \ is a known loading boundary 
which could be written using physical Cauchy stress. 



4.1 Finite element approximation 

Through this paper, we use a three dimensional finite clement with Lagrange trilincar interpolation for 
the displacements and uniform Lagrangian multipliers to compute an approximate solution of Eqs.(14)- 
(15)-(16) on a rectangular mesh (see Fig. 4.1), where we neglect the acceleration and body forces (b = 0, 
7 = 0). This element is commonly used and is relevant for the finite clement approximation of this type 
of problem where kinematics constraints must be satisfied (for more details see [18, 9, 10, 11, 5, 22]). 




Position 
#Aj #a) 



1 Y, 



Figure 4: Ql — QO element for displacements, pressure and pseudo-active Lagrange multiplyer. 

Let (£k) the Lagrangian normalized finite element coordinates (Figure A), the deformed geometric coor- 
dinates 9 a in clement e are interpolated as 



0"= E VVe)(£l,6,6)C(e) 
n(e) = l 



(18) 



where Vn(e) is the base function associated with the local node n(e) and 0",s is the a-coordinate of the 
local node n of element e. 



Let fi^ 6 -* be the connectivity matrix defined by 



n n(e) = f 1 if A(n(e),e) = A 
A 10 otherwise 



(19) 



The FE approximation of Eqs.(14)-(15)-(16) is 



E E °A 6) / ^ J *> a V/(Vn(e))d^ = (l-^H)E E "a 6) / ^ ^n(e) ^ 



(20) 



e n(e) = l 



e rt(e)=l 



Ve , y (det 5 $ - l) dV = 



(21) 



/ 



Ve , 5 ah / /i(/ 4 ,/ 6 )^ = 



(22) 



with A = 1, • • • ,A max , a — 1,2,3, where V e is the volume of the element e, A?. e is the part of A e 
(boundary of the element e) non subject to displacement conditions. 



4.2 Finite element solution method 

We proceed in two steps. The first one consists in the determination of the pseudo- active stresses T{ 
and T°£ as functions of the activation parameter (3 G [0,1] by looking for the state A (Sah = 1)- We 
solve the system (20)-(21)-(22) with zero right hand side for (20) (free active contraction) and P 1 / given 
by Eq.(26) . The unknowns of this nonlinear system of equations are (#^,pA(e), (kt(e)) with a — 1,2,3, 
A = 1, • • • , A. max and e = 1, . . . , e max where e max is the total number of elements involved in the mesh. 
We derive T' A and T c / for a given (3 according to Eq.(12). Then in a next step we can compute any 
physiological active loaded state C solving the system (20)-(21) with Sah = and P^f given by Eq.(27). 
To solve the sytem in both cases we use the Powell method [21] implemented in the package MINPACK 
[17]. 
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5 Results and discussion 



This section is devoted to the numerical simulation of two types of material: 

• a thin sample of living myocardium for which a cartesian coordinate is used, 

• an active thick-walled cylinder for which cylindrical coordinate is used. 

For these two simple configurations the exact displacements are solutions of a nonlinear system and can 
be computed with a high degree of accuracy. A very good agreement between the exact and the computed 
solutions of the finite element nonlinear system (20)-(21)-(22) is obtained. More precisely in all the cases 
the L 2 norm of the error is less than 1CP 09 . 

5.1 Case of a thin sample of living myocardium 

We simulated the loading of a thin sample of living myocardium (1.0 x 1.0 x 0.1 cm 3 ) in which the MF 
are uniformly oriented in one direction. The coefficients involved in the strain energy-function are those 
of Lin and Yin [15]: Cf =0.292 kPa, Cf =0.321, Cf =-0.260, C\ =0.201, C?=-3.870 kPa, Cf =4.830 kPa, 
Cf =2.512 kPa and C.4=0.951 kPa. For the beating tension, a good agreement between the previous 
experimental results and our theoretical solution is obtained for T^=0.6 kPa. Nevertheless, the control 
simulation was performed with a/D=0.2 and T(°)=35 kPa. This higher value of is more adapted to 
the description of the left ventricular performance [19]. 

Influence of the collagen network on the systolic wall thickening- The free contraction test is 
performed with no external displacement or force on the boundaries of the sample, but just in activating 
the tissue. In this simulation we used the following activation function : /3(s) = sin 2 (its). Compare to the 
case where the kinematic constraint is not taken into account, one can see an increase of the cross-fiber 
extension ratio which is in the tangential plane of the ventricular wall (Figure 5). At the end-systolic 
state (i.e. when (5 — 1), this ratio goes from the value 1.25 if we neglect the coupling effect, to 1.45 when 
considering the kinematic constraint induced by the collagen. So, the connective tissue could account 
for 16 % of normal end-systolic wall thickness. This increase is clearly dependent of the geometrical 
parameter ratio a/D and the maximal beating tension T( ). 
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Figure 5: Free contraction test with (3(s) 



2 (7rs): effect of the pseudo-active kinematic constraint. 



The empty and full symbols indicate that the coupling effect is acting or not, respectively. The fiber and 
cross-fiber directions are noted (f), (cf) and (cf ') and are defined in figure 2. Arrows show the curve 
modification when the pseudo-active kinematic constraint behaves. 



Influence of the collagen network on the pseudo-active tension- Table 1 shows the effect of the 
geometrical parameter a/D and the maximal active tension T^-°\ on the fiber and cross-fiber stresses 
(noted (Jn and tr 2 2, respectively). These effects were given in the case of an equibiaxial extension loading 
(Xf = X c f = 1.2) of an activated sample of myocardium (f3 = 1). These two stresses increase with 
T", but are not very sensitive to the geometrical ratio a/D. We can observe also, that by neglecting 



the interaction between the collagen network and the MF: (i) the cross-fiber stress is not affected by 
the amplitude of the beating tension, and (ii) the stress ratio 022/0-11 decreases when increases. 
These results mean that the usual strain-energy functions considered for the myocardium are not able 
to generate any transverse pseudo-active tension. Moreover, the results obtained for the uniaxial tests of 
an active or a passive sample, with or without the effect of the collagen on the MF, are shown in Figure 
6. Because the coupling effect between the collagen and the MF is an active mechanism, the passive 
stress-strain relations are not affected by the kinematic constraint. The mechanical properties of the 
active tissue, in the fiber and cross-fiber directions, become comparable when the coupling effect acts. 



Table 1: Effect of active tension > and geometrical parameter a/D 



a/d 




5 


15 


25 


35 


45 


0.10 


0-22/0-11 (%) 


50.92 


53.73 


57.65 


60.25 


61.90 


022 (KPa) 


5.94 


11.44 


17.78 


24.32 


30.88 


0.15 


0-22/0-11 (%) 


51.44 


54.56 


58.60 


61.30 


63.01 


0-22 (KPa) 


5.97 


11.48 


17.80 


24.31 


30.84 


0.20 


0-22/0-11 (%) 


51.94 


55.33 


59.47 
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64.04 


0-22 (KPa) 


6.00 


11.51 


17.79 


24.25 


30.73 


No kinematic constraint 


0-22/0-11 (%) 


36.78 


19.88 


13.62 


10.36 


8.36 


0-22 (KPa) 


4.32 


4.32 


4.32 


4.32 


4.32 




Figure 6: Active and passive uniaxial extension tests: effect of the pseudo-active kinematic constraint. 
The empty and full symbols indicate that the coupling effect is acting or not, respectively. The fiber and 
cross-fiber directions noted (f) and (cf) are defined in figure 2. Arrows show the curve modification when 
the pseudo-active kinematic constraint behaves. 



5.2 Case of an active thick-walled cylinder 



We simulate the mechanical behaviour of an active artery under physiological blood pressure Pint ■ This 
artery is modelled by a thick-walled cylinder with internal radius Ri n t — 2 mm, external radius R ex t = 
3.5 mm and height L — 2 cm. We assume that the medium is made of a hyperelastic anisotropic material 
with fibers oriented in the circumferential direction. This simulation does not take into account all the 
complexity of the structure of an artery, so the following results must be viewed as a first approach and 
we focus on some qualitative aspects, particularly the wall thickening effect. 




Figure 7: a simplified model for the spatial organization of the collagen fibrils. The fibrils are oriented 
in the z-direction (c/) and in the right part of the figure we suggest that the fibers are oriented in the 
circumferential direction 

In 1902, Bayliss suggested that the distension of the vessel by blood pressure could act as a mechanical 
stimulus to the vascular smooth muscle cells, thereby contributing to their tone [1]. However, conclusive 
experimental support for this concept was available only recently. We now know that the degree of vascu- 
lar distension appears to be a factor of importance in determining vascular tone. We used the suggested 
constitutive law to model a hypothetical autoregulation mechanism. 

For this simulation, the active fiber tension as well as the rheological change are in phase with the pul- 
satile pressure, and we use as input data the following functions: (3(s) = sin 2 ns with Pint = 8 + 10 sin 2 irt 
(Fig. 8). The resulting variations of the thick- walled cylinder radii are presented in Figures 10-11 for this 
autoregulation law based on fluid pressure. The autoregulation is defined as the relationship between the 
activation function (3(s) and the pulsatile blood pressure Pint- Very interestingly, the results show that 
the kinematics of the arterial wall may be more sensitive to the change of mechanical properties than to 
the blood pressure. In other words, it appears that if there is no kinematic constraint due to the fibrils, 
the internal and external radii increase when the blood pressure decreases. In fact, during this decrease 
of pressure, we assume that the material becomes more compliant. Thus, the wall kinematics is mainly 
driven by the change of rheology. Furthermore, although the pressure and activation are in phase, you 
can create with this autoregulation law some delay in the kinematic response. Therefore we believe that 
the pressure-activation interaction is a fundamental mechanism which must be well modeled to describe 
accurately the behaviour of the arterial wall under physiological or pathological conditions [28, 24]. 
On the other hand, when the kinematic constraint is activated, we observe that the interior radius still 
decreases but less than previously, due to the residual constraints T A and (see Fig. 10). The more im- 
portant effect is observed for the exterior radius (Fig. 11), due to a greater residual constraints T;{ and T c / 
and the absence of exterior stress. The global result is a wall thickening effect and a contraction in the 
z-direction due to the incompressibility condition (see Fig. 9) 
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Figure 8: Pressure and tonus as a function of the degree of activation 
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Figure 9: height of the cylinder as a function of the degree of activation 
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Figure 11: exterior radius of the cylinder as a function of the degree of activation 



6 Conclusion 

This study shows that the connective tissue skeleton in the normal and pathological left ventricle may 
have a large influence on the cardiac performance. A new constitutive law has been developed for large 
deformations of an incompressible hyperelastic, and anisotropic living myocardium. This work is based 
on the idea that the connective tissue is physically coupled to the muscle fibers which seems reasonable 
with regard to the available observations. Nevertheless, additional experimental works must be done in 
order to support this assumption and to study thoroughly the spatial organization of the myocardial 
collagen fibrils under normal and pathological conditions. 



A Coordinate systems 
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Figure 12: Coordinate systems (adapted from Costa et al. [6]). 



Table 2: Notations for the coordinate systems used to formulate the finite element method (adapted from 
Costa et al. [6]). (I) Rectangular cartesian reference coordinates, (II) Curvilinear world coordinates, (III) 
Normalized finite element coordinate, (IV) Locally orthonormal body/fiber coordinates (adapted from 
Costa et al. [6]). 
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B Second Piola-Kirchoff stress tensor at states A and C 

In the case of the free contraction state A, using Eqs.(4)-(13)-(10)-(17), we can write the components 
P^ J of the second Piola-Kirchoff stress tensor in the state A base tensor (G^ ® G^)a under the 
form: 



P\ J = - P A9 (X)IJ + 2G^ IJ Wt + {2Wl-q A h i ).fP I .f { ; )J 

while in the case of the active loaded state C, using Eq.(9)-(17), the components Pq J of Pc in the state 
C base tensor (G^ ® Gj')c are 



P£ 7 = -p c g (x)IJ + 2G^ ij W^ + 2WUp ]I !p 



8W* F)W f)W f dW pseudo 

where W* = ^- = ^fi + /? ^ + S AH i = 1, 4 (23) 

and = i = 4) 6 (24) 

/jf^ 7 , fp^ 1 are respectively the components of the unit vectors fp, fp- in the base (Gj X \ / = 1,2,3) 
and /jy ^ 7 , f^ 1 with il = A or C, are respectively the components of the unit vectors f#, i' H in the base 



(s\ X \ I = 1) 2, 3). The metric tensors G^ IJ , g^ IJ are defined in table 1. 

Following the definition of the locally orthonormal body/fiber coordinate system we have /p X ^ J = 5 11 and 
fp^ 1 = S 21 . On the other hand the vectors f# and are respectively defined through: 

tff HWpII || fj?>' g <*> |l IIW^II Hf^g^ll 

thus = pr — and = j—, — and we get finaly: 

II Si II " II 82 II 

Pj J = -p J 4 5 (a:)/J + 2G( :i; ) /J T^ 1 * + (2Ty 4 *- gj 4/i4)<5 1 ^ 1J 

11 sr 11 2 



and 



P C J = -p c ff ^ IJ + 2G {x)IJ W* + 2 Wl 5 U S 1J 

+(/?T<°> + Tj) + I*' (27) 

II II 2 II g^ II 2 
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